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Using high-resolution particle image velocimetry we measure velocity profiles, the wind 
Reynolds number and characteristics of turbulent plumes in Taylor-Couette flow for a 
radius ratio of 0.5 and Taylor number of up to 6.2 • 10^. The extracted angular velocity 
profiles follow a log-law more closely than the azimuthal velocity profiles due to the 
strong curvature of this = 0.5 setup. The scaling of the wind Reynolds number with 
the Taylor number agrees with the theoretically predicted 3/7-scaling for the classical 
turbulent regime, which is much more pronounced than for the well-explored rj = 0.71 
case, for which the ultimate regime sets in at much lower Ta. By measuring at varying 
axial positions, roll structures are found for counter-rotation while no clear coherent 
structures are seen for pure inner cylinder rotation. In addition, turbulent plumes coming 
from the inner and outer cylinder are investigated. For pure inner cylinder rotation, the 
plumes in the radial velocity move away from the inner cylinder, while the plumes in 
the azimuthal velocity mainly move away from the outer cylinder. For counter-rotation, 
the mean radial flow in the roll structures strongly affects the direction and intensity of 
the turbulent plumes. Furthermore, it is experimentally confirmed that in regions where 
plumes are emitted, boundary layer profiles with a logarithmic signature are created. 

Key words: 


1. Introduction 

The paradigmatic Taylor-Couette (TC) flow has long been a flow configuration of great 
interest to fluid dynamicists. It consists of flow between two coaxial cylinders that can 
independently rotate, see figure 1. This system has been used extensively as a model in 
fluid physics because it is a closed system, has a relatively simple geometry and therefore 
has multiple symmetries. After early investigations (Couette 1890; Mallock 1896; Taylor 
1923), Wendt (1933) started studying turbulence in this system. For a further historical 
overview, the reader is referred to Donnelly (1991). For a review on TC flow at the 
onset of instabilities and slightly above, see Fardin et al. (2014). The state-of-the-art of 
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Figure 1. (a) Schematic of the vertical cross-section of the Taylor-Couette setup, showing the 
horizontal laser sheet which is imaged from above trough the transparent top plate. The laser 
can be traversed vertically, (b) Horizontal cross-section of the setup. The red rectangle represents 
the typical field of view. 


high Reynolds number Taylor-Couette turbulence is treated by Grossmann, Lohse & Sun 
(2016). 

The two geometrical control parameters of TC flow are the ratio 77 of the inner and 
outer cylinder radius and the aspect ratio F = djL^ where d is the gap width and 
L the height of the cylinders. In this work we use a small radius ratio of 77 = 0.5, 
corresponding to a relatively wide gap. The vast majority of existing work (see overview 
by Dubrulle et al (2005) and review by Grossmann et al (2016)) focuses on a radius 
ratio of 0.71 or higher. Lower 77 experiments concern the Rayleigh-stable regime (Ji et al. 
2006), mean flow and turbulence characteristics (van Hout & Katz 2011) and global 
torque measurements (Merbold et al. 2013). In addition, numerical work has been done 
for 77 = 0.5 (Dong 2007; Ghouippe et al. 2014; Ostilla-Monico et al. 2014c). The radius 
ratio is a key control parameter (Eckhardt et al. 2007) in TG flow and it has been found 
that a low radius ratio Taylor-Gouette system behaves differently as compared to higher 
77 setups in several aspects, which will be outlined below. The underlying reason for 
the different behaviour is the strong boundary layer asymmetry. Because the ratio of 
inner and outer boundary layer thicknesses scales as 77 ^ (Eckhardt et al. 2007), a strong 
asymmetry between the inner and outer boundary layer exists for small radius ratio. 

As dimensionless control parameters of the system we use the (negative) rotation ratio 
a = —ujol^i and the Taylor number (Eckhardt et al. 2007): 


_ (1 + rj^ {tq - Tirin + ro)‘^{L0i - uior 

64772 jy 2 


where uji^o are the angular velocities of the inner and outer cylinder, Vi^o are the inner 
and outer cylinder radii, 77 = r^/ro is the radius ratio, and u is the kinematic viscosity. 
The driving of the system can also be described by the shear Reynolds number (Dubrulle 
et al. 2005): 
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When increasing the driving strength (i.e. the Taylor number), TC flow first gradually 
undergoes a transition from a purely azimuthal, laminar state to a state where the bulk 
becomes turbulent while the boundary layers still remain laminar. The latter state is 
called the classical regime (Grossmann & Lohse 2000). By further increasing Ta, the 
ultimate regime (Kraichnan 1962; Grossmann & Lohse 2011; He et al. 2012 ) is reached, in 
which also the boundary layers are turbulent. The signature of these turbulent boundary 
layers are logarithmic velocity profiles which have recently been found for a radius ratio 
of 77 = 0.716 (Huisman et al. 2013) and for 77 = 0.909 (Ostilla Monico et al. 2015). 
Grossmann et al. (2014) have derived theoretical velocity profiles, and found that the 
angular velocity profile follows a log-law more closely than the azimuthal velocity profile, 
an effect that is more pronounced with the stronger asymmetry for smaller 77 . In section 
3.1 we set out to investigate the correspondence of experimentally measured angular and 
azimuthal velocity profiles to a log-law for 77 = 0.5 and Taylor numbers up to 6.2 • 10^ 
at the onset of the ultimate turbulent regime. Furthermore, the velocity profiles are 
compared to existing experimental and numerical work. 

In addition to affecting the velocity profiles, the radius ratio strongly influences the 
transitional Taylor number for the ultimate regime (Ravelet et al. 2010; Merbold et al. 
2013; Ostilla-Monico et al. 2014a,c). For 77 = 0.5 the ultimate regime does not start before 
Ta = 10^^ (Merbold et al. 2013; Ostilla-Monico et al. 2014c), whereas the transition for 
a higher radius ratio of 77 = 0.71 already occurs at Ta = 5 • 10^ (van Gils et al. 2012 ; 
Ostilla-Monico et al. 2014 c). Very different scaling of the angular momentum transfer 
and the “wind” in the gap with the driving paramater Ta exist for the classical and 
ultimate regime (Grossmann & Lohse 2000, 2011). The late onset of the ultimate regime 
for 77 = 0.5 gives us the opportunity to analyze the scaling of the wind Reynolds number 
with Taylor numbers in the classical TG regime, up to the onset of the ultimate regime 
and compare it with theoretical predictions (section 3.2). 

Another area where the influence of the radius ratio is apparent, is the angular 
momentum transport trough the gap. This transport is a key response parameter of the 
system and is directly connected to the torque required to maintain constant cylinder 
velocities (Eckhardt et al. 2007). The radius ratio strongly influences the rotation ratio 
for which optimal momentum transport occurs (van Gils et al. 2012 ; Merbold et al. 2013; 
Ostilla-Monico et al. 2014a). For example, the optimal momentum transport occurs at 
<^opt ~ 0.33 for T] = 0.714 (van Gils et al. 2012 ) whereas it is Uopt ~ 0.20 for 77 = 0.5 
(Merbold et al. 2013). 

Roll structures play an important role in the momentum transport in the Taylor- 
Gouette system (Fenstermacher et al. 1979; Andereck et al. 1986; Lewis & Swinney 1999; 
Huisman et al. 2014; Ostilla-Monico et al. 20146). With increasing Taylor number, the 
system undergoes transitions from a purely laminar state to having steady, coherent 
Taylor vortices, to having modulated, wavy Taylor vortices which eventually evolve into 
chaotic turbulent Taylor vortices (Lewis & Swinney 1999). But even in the ultimate 
regime, for some a > 0 , it was found that remnants of these rolls are present in time- 
averaged quantities (Tokgoz et al. 2011; Huisman et al. 2014; Ostilla-Monico et al. 20146), 
whereas for a = 0 (pure inner cylinder rotation) these rolls vanish for large Taylor number 
(Lathrop et al. 1992). In section 3.3 we will characterise the roll structures for 77 = 0.5 
at both pure inner cylinder rotation a = 0 and slight counter-rotation a = 0 . 2 , at which 
optimal transport occurs. 

It was recently proposed that the logarithmic velocity profiles in ultimate TG flow 
are triggered by turbulent plume ejection (Ahlers et al. 2014; Ostilla-Monico et al. 
20146; van der Poel et al. 2015), which in turn is regulated by the relative strength 
of the axial and radial mean flow. Using time-resolved azimuthal and radial velocity field 
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Figure 2. The explored parameter space. The Taylor number (Ta) is varied for inner cylinder 
rotation (a = 0), and axial scans have been performed for (negative) rotation ratios a = 0 and 
a = 0.2. 

measurements at several heights in the setup we will investigate the connection between 
the roll structures that generate strong radial and axial flow, and the turbulent plumes 
that emanate from either the inner or outer cylinder (section 3.4). Futhermore, in section 
3.5 we will quantify how these plumes affect the logarithmic nature of the velocity profiles. 

2. Setup & explored parameter space 

2.1. Setup 

The experiments were carried out in the Cottbus Taylor-Couette facility (Merbold 
et al. 2011, 2013). The inner and outer cylinder radii are = 35.0 ± 0.2mm and = 
75.0 ± 0.2mm, respectively, the height of the setup is L = 700 mm. This gives a radius 
ratio of r] = 0.5 and an aspect ratio of T = 20. The cylindricities of the cylinders that 
were used are 0.4 mm and 0.3 mm for the inner and outer cylinder, respectively. The 
maximum rotation rates are 5 Hz for both the inner and outer cylinder. 

The end plates rotate with the outer cylinder. The outer cylinder and the top plate 
are transparent, making the setup ideally suited to be used in combination with particle 
image velocimetry (PIV). A high-resolution PIV camera (LaVision Imager sCMOS) 
with a resolution of 2560 x 2160 pixels and a frame rate of 50 Hz is installed above 
the top end plate, pointing downwards. Water is used at the working fluid (20°C, u = 
1.0 X 10“^m^/s). The water contains fluorescent particles (Dantec Dynamics, PMMA- 
RhB, 1-20 jam) with a maximum Stokes number of St = r^/r^ ~ 10“^ <C 1, which 
means that they faithfully follow the flow and can be considered as tracer particles. 
The flow is illuminated by a horizontal light sheet from a high-powered pulsed NdiYLF 
dual cavity laser (Litron LDY303HE). Because of the high-resolution PIV camera, very 
high resolution measurements of the flow fields can be achieved. The imaging of the 
full width of the gap combined with a vector grid of 16 x 16 pixels with 50% overlap 
results in a velocity vector spacing of 0.13 mm. The PIV system is operated in dual frame 
mode, allowing for an interframe time At smaller than the inverse frame rate. The PIV 
image pairs are processed using LaVision DaVis software, after which the flow fields are 
transformed to the radial velocity Ur{0, r, 2 ;, t) and the azimuthal velocity ue{0, r, z, t). 

2.2. Explored parameter spaee 

The parameter space that was explored for this work is shown in figure 2. The first set 
of measurements that is treated here consists of varying the Taylor number between Ta = 
5.8 • 10^ and Ta = 6.2 • 10^ (corresponding to a shear Reynolds number of Res = 6.0 • 10^ 
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Figure 3. Overview of the flow profiles for varying Taylor number with pure inner cylinder 
rotation a = 0 at mid-height h/L — 0.5. (a) Snapshot of the flow field for Ta = 6.2 • 10^. 
The colours and lengths of the arrows indicate the norm of the velocity |u|. (b) Averaged 
flow field over 10000 PIV image pairs at Ta = 6.2 • 10^. (c) Radial velocity profiles across the 
gap of the TC apparatus, normalised by the inner cylinder velocity. For lower Taylor numbers 
there is still a strong radial flow, which can be attributed to the presence of Taylor vortices, 
(d) Azimuthal velocity profiles normalised by the inner cylinder velocity. For increasing Taylor 
number, the profiles become flatter and the boundary layers steeper, (e) and (f) Magnification 
of the azimuthal velocity profiles in (d) close to the cylinders. 


to 6.2 • 10^) for inner cylinder rotation only (a = 0) at mid-height {h = 0.51/). These 
measurements were performed as a single continuous experiment; the cylinder velocity 
was increased slowly between experiments and before each measurement approximately 
10 minutes were taken to give the flow ample time to stabilise. Each measurement in this 
set consists of 10000 PIV image pairs, which were recorded at 25 frames per second. This 
corresponds to 235 and 2418 cylinder revolutions for Ta = 5.8 • 10^ and Ta = 6.2 • 10^, 
respectively. 

In addition to investigating the Taylor number dependence, the height of the laser 
sheet was varied to characterise the axial dependence of the flow field. At 13 different 
heights with 7.5 mm spacing, 5000 image pairs of the velocity field were recorded. This 
was done at Ta = 4.2 • 10^ for both pure inner cylinder rotation and optimal (Merbold 
et al. 2013) counter-rotation a = 0.2. 


3. Results 

3.1. Azimuthal and angular velocity profiles 

In this part the properties of flow profiles for varying Taylor number will be inves¬ 
tigated. An example of an instantaneous flow field is shown in figure 3(a), while the 
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Huisman et al. (2013), 77 = 0.716: 

Chouippe et al. (2014), 77 = 0.5: 

Ta = 5.8*10^ 

Ta= l.MO^ 



Ta = 9.9 10 
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Ta = 2.M0^ 
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Ostilla-Monico et al. (2014c), 77 = 0.5: 

Ta = 4.4-10^ 
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Ostilla-Monico et al. (2015), 77 = 0.5: 
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Ta = 6.210^^ 
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Figure 4. Velocity profiles for inner cylinder rotation a = 0 and varying Taylor number, (a) 
Azimuthal velocity profiles near the inner cylinder (f G [ 0 , 1 / 2 ]). The velocity u'^ is defined as 
{u 0 {ri) — U 0 {r))/ur,i with the friction velocity Ur,i = \^rw,ilp containing the wall shear stress 
Tw,i = Tj^'Kr^L with the torque T. The wall normal distance is scaled using the viscous length 
scale — v (b) Azimuthal velocity profiles near the outer cylinder (f G [ 1 / 2 , 1 ]). For the 
outer cylinder case, the velocity vV is defined as {u 0 {r) — U 0 {ro))/ ut,o with Ut,o = \Jtw,oI p and 
Tw,o — Tl2'Kr‘lL. The normalised wall normal distance is — {tq — r^jd^.o with — njuT^o- 
(c) Angular velocity profiles near the inner cylinder (f G [0,1/2]). The velocity is defined as 
{uj{ri) — uj{r))/{uT,i/ri) (d) Angular velocity profiles near the outer cylinder (f G [1/2,1]). For 
the outer cylinder case, the velocity is defined as (ct;(r) — uo{ro))/{ ut^o/to) All figures include 
the logarithmic law of the wall from Prandtl and von Karman = 1/nhiy^ B with the 
typical values oi n — 0.40 and B — ^.2 (see Marusic et al. (2013) and references therein), the 
viscous sublayer u'^ = DNS data from Chouippe et al. (2014), Ostilla-Monico et al. (2014c) 
and Ostilla Monico et al. (2015) at 77 = 0.5 and measurement data from Huisman et al. (2013) 
at ?7 = 0.716. 
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averaged field is shown in (b). From (c) it is apparent that at this specific axial position 
(mid-height) a significant mean radial flow still exists for the lower Taylor numbers up to 
approximately 10^. These are Taylor vortices, which at a = 0 disappear for higher Taylor 
numbers (Lathrop et al 1992). With increasing Taylor number the azimuthal velocity 
profile in the bulk becomes flatter (figure 3(d-f)), but it is still apparent that the strong 
curvature of a radius ratio of 77 = 0.5 creates a significant asymmetry between the inner 
and outer boundary layer. 

In figure 4 profiles up to the middle of the gap are shown, normalised to the wall normal 
distance from both the inner and outer cylinder, and to the azimuthal velocity u'^ 
((a) and (b)) and angular velocity co’+ ((c) and (d)). The friction velocity Ur that is used 
in this normalisation is calculated with global torque data from Merbold et al (2013). 

As the driving is increased, the profiles slowly approach the Prandtl-von Karman log- 
law, although at these Taylor numbers the log layer is not yet fully developed. For 
comparison, profiles at higher Taylor numbers and aspect ratio of 77 = 0.716 from 
Huisman et al. (2013) are also plotted. The data show good agreement with direct 
numerical simulations (DNS) at the same radius ratio 77 = 0.5 from Chouippe et al. 
(2014) at Ta = 1.8-10^ Ostilla-Monico et al (2014c) at Ta = 1.0-10^ and Ta = 1.0-10^^ 
and Ostilla Monico et al (2015) at Ta = 1.1 • 10^^. 

When comparing the azimuthal velocity 77 + to the angular velocity (i.e. (a) to (c) 
and (b) to (d)) it can be seen that the angular velocity profiles curve upwards more 
compared to the azimuthal ones. This is consistent with the theoretical argument of 
Grossmann et al. (2014) that the angular velocity profile co’+ is closer to a log-law than 
the azimuthal velocity profile 77 +. The effect is much more pronounced for the smaller 
77 = 0.5 here as compared to the 77 = 0.716 from Huisman et al. (2013). Note that, 
for example, the lower Taylor number 77 + profiles at 77 = 0.716 are above the profiles 
measured here, but that they cross when represented as This clearly shows the 
influence of the large curvature for small radius ratio TC setups. 

3.2. Wind Reynolds number 

The degree of turbulence of the wind in the gap of the cylinders, which measures the 
strength of the secondary flows Ur and 77 ^, can be characterised by the wind Reynolds 
number. We use the standard deviation of the radial velocity to quantify the wind 
Reynolds number: 

RCw = G{Ur)d/v (3.1) 

with the gap width d = To — r^. In the analogy between TC and Rayleigh-Benard flow, the 
unifying theory of Grossmann & Lohse (2000) predicts that the wind Reynolds number 
in the classical turbulent regime scales as: 

Rew oc Ta^/^. (3.2) 

In contrast, in the ultimate turbulent regime the scaling was found to be Re^ oc Ta^'^^ 
(Grossmann & Lohse 2011), which was confirmed experimentally by Huisman et al. ( 2012 ) 
for TC flow with 77 = 0.716 for Ta up to 6.2 • 10^^. 

The highest Taylor number we achieve in the present measurements is Ta = 6.2 • 10^, 
which for 77 = 0.5 with the higher transitional Taylor number of Ta = 10^^ (Merbold et al 
2013; Ostilla-Monico et al. 2014c) implies that we are in the classical turbulent regime. 
Therefore we expect the theoretically predicted 3/7-scahng to hold. This classical scaling 
was measured in Rayleigh-Benard convection (He et al. 2012 ), but to our knowledge has 
not yet been confirmed for TC flow. 

The wind Reynolds number Re^ is extracted from the velocity field 77 ^(d, r, z, t) by 
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Figure 5. The scaling of the wind Reynolds number with Taylor number for inner cylinder 
rotation a = 0. (a) The standard deviation of the radial velocity over 10000 frames (400 s) 
and the azimuthal direction, normalised by the inner cylinder velocity, as a function of the 
normalised radial position. See the legend in figure 3 or 4 for the values of the Taylor numbers. 
The three areas A (f G [1/8, 3/8]) , B (f G [3/8, 5/8]) and C (f G [5/8, 7/8]) indicate regions 
over which a is averaged, (b) Rcw versus Ta averaged over f-values corresponding to regions 

A, B, C and all three combined (r G [1/8, 7/8]). In addition, data from Huisman et al. (2012) 
with their fit are shown. This data uses the equivalent to area B for averaging, (c) Data in (b) 
compensated by Ta^^^. The fits are Rcw = 0.162Ta°‘^^^ for area A, Rcw = 0.193Ta°'^^^ for area 

B, Rcw = 0.240Ta°‘^°^ for area C and Rcw = 0.192Ta°‘^^^ for the three areas combined. 


calculating the standard deviation a{ur{r,0A)) over time and the azimuthal direction, 
which is then averaged over a certain range in the gap, see figure 5. When looking 
at the shape of the r-dependence of a, it is apparent that an asymmetry is present. 
The fluctuations exhibit a maximum around a quarter of the gap width away from the 
inner cylinder. It is likely that the strong curvature of the 77 = 0.5 setup causes this 
assymmetry. Calculating the wind Reynolds number thus comes with an arbitrariness 
in the area to choose for averaging the profile of cr(r), although excluding the boundary 
layers is reasonable because the interest lies in the convective bulk transport, as opposed 
to the diffusive transport in the boundary layers. The dependence of Rew(Ta) is shown 
in figure 5(b) for different regions of averaging. 

From figure 5(c) it can be seen that a clear power law scaling is present with an 
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Figure 6. Dependence of flow profiles on axial position for rotation ratios a = 0 (a-d) and 
a — 0.2 (e-h) for Ta = 4.2 • 10^. The normalised radial velocity is Ur — Ur(r)/{ui — Uo) and 
the azimuthal velocity is uq — {ue(r) — Uo)/{ui — Uo). The data is represented in a colour map 
(a,c,e,g) and as prohles in (b,d,f,h). For the visualisation of the colour map, bilinear interpolation 
is used. 

exponent that ranges from 0.402 to 0.434, depending on the region of averaging. This 
is consistent with the theoretically predicted classical regime scaling of 3/7 ^ 0.429 
(Grossmann & Lohse 2000). 
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3.3. Roll structures 

The phenomena of roll structures, turbulent plumes and logarithmic velocity profiles 
in Taylor-Couette flow are intimately connected (Ostilla-Monico et al. 20146,c). We will 
first show an example of roll structures that are observed in TC flow. Flow profiles at 
several heights for both a = 0 and a = 0.2 are measured, allowing us to visualise roll 
structures as shown in figure 6. At this reasonably high Taylor number of Ta = 4.2 • 10^, 
no clear roll structures exist for inner cylinder rotation only. The radial velocity shows 
some very faint indications of rolls, with a maximum value of Ur = 0.01 in the center of 
the gap. In contrast, in the case of counter-rotation (a = 0.2), very strong roll structures 
can be seen. The fact that there are strong rolls at optimal counter-rotation and no 
clear structures at inner cylinder rotation only, is corroborated by recent work for higher 
radius ratio (Tokgoz et al. 2011; Huisman et al. 2014; Ostilla-Monico et al. 20146). 

In the radial flow, between two rolls at approximately hjL = 0.51, there is very strong 
outwards radial flow (see figure 6(e)). Specifically, the height 0.52L corresponds to the 
bottom of a roll, while 0.50L corresponds to the top of the roll below. Exactly in between 
the rolls at 0.511/, the radial flow surprisingly exhibits a local minimum. In figure 6(g,h) 
it can be seen that in the case of positive (outward) radial flow, inner cylinder velocity 
is advected outwards, shifting the profiles upwards at the inner half of the gap. For 
negative (inward) radial flow, the profiles shift downwards in the outer half of the gap. 
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In summary, the roll structures are represented by a large secondary flow of Ur^ which 
advects velocity from the cylinders and changes the azimuthal velocity profiles. In the 
following two subsections we will elaborate on the characteristics of turbulent plumes 
and their effect on the velocity profiles. 

3.4. Turbulent plumes 

In turbulent Taylor-Couette flow, structures detach from the boundary layers, called 
herring-bone streaks (Dong 2007) or velocity plumes (Ostilla-Monico et al 20146), 
in correspondence to the thermal plumes in the analogous turbulent Rayleigh-Benard 
convection. These plumes in TC flow are large-scale spatial and temporal fluctuations in 
the velocity fields that detach from either the inner or outer cylinder and advect velocity 
from the respective cylinder. Because of the sufficiently high frame rate of 25 fps and 
large amount of data recorded in our experiments, these plumes can be resolved in time, 
and the typical velocity can be extracted. We are interested in understanding how these 
plumes are affected by roll structures and whether some scaling relation exists for the 
typical plume velocity as a function of the driving parameter Ta. 

The time dependence of the profiles of the azimuthal and radial velocity can be 
represented by the spatio-temporal fields Ur{r,t) and uo{rA), with the value of the 
velocities represented by a colour map, as shown in e.g. figure 7(a). In this way, structures 
or fluctuations of the profiles that propagate inwards or outwards become visible as 
diagonal lines. It has to be stressed that these diagonals do not represent a nonzero 
mean radial flow, but only the fluctuations on these profiles that travel outward or 
inward. In order to assign a velocity to these plumes we employ an analysis which is 
based on applying the following afflne shearing transformation to the functions i4^(r, t) 
and uo{r, t): 

(t, r) (t + d • r/vb, r) (3.3) 

with Vb the plume velocity. For a certain amount of shearing, the plumes which were 
represented by diagonal lines, become vertical lines. To And this optimal value of shearing, 
we postulate that when averaging the sheared function over f, the standard deviation 
of the resulting time signal has a maximum as a function of the shearing value. At 
this optimal shearing, the plumes are represented by narrow peaks in the time signal, 
creating a high standard deviation. Using a golden section search algorithm (Kiefer 1953) 
the optimal shearing value is found, which corresponds to the mean velocity of the plumes 

Vb. 

The dependence of the plume velocity on the axial position and the rotation ratio is 
shown in figure 7. Again, for inner cylinder rotation only (a = 0), no height dependence 
can be seen. The fluctuations on the radial velocity are moving outwards with a velocity 
'^6/('^i —Uo) = 0.25 ±0.02 {uo = 0, 2(7 error). When the azimuthal velocity is considered, 
the picture looks different. Crossing diagonal lines can be seen, meaning that plumes 
emit from both the inner and outer cylinder. Despite of this, a clear peak was found 
for inward plumes with an average relative velocity of Vbl{ui — Uq) = 0.20 ± 0.02. The 
velocity of these plumes is of the same order, but it remains to be explained why the 
main plume direction is different for the two velocity components. 

For counter-rotation (a = 0.2), the bottom half of figure 7 clearly indicates that roll 
structures are present. There is a correspondence between the sign of the mean radial 
velocity and the direction of these turbulent plumes, coming from either the inner or 
outer cylinder. Compare figure 6(e) with figure 7(h,k). For both velocity components, 
the direction of the plumes corresponds to the direction of the mean radial velocity. The 
mean radial velocity forces the plumes strongly in the inward or outward direction, also 
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Figure 7. Overview of turbulent burst velocity depending on axial position, for both pure inner 
cylinder rotation a = 0 (a-f) and counter-rotation of a = 0.2 (g-1), at Ta = 4.2 • 10®. The radial 
velocity dependence is shown in (a,b,c,g,h,i) and the azimuthal velocity in (d,e,f,j,k,l). (a,d,g,j) 
Example of the radial and azimuthal velocity as a function of the normalised gap distance and 
time for the axial position of 0.50L. Only 10 s of the total 200 s per experiment are shown for 
clarity. Diagonal lines with a positive (negative) slope correspond to plumes coming from the 
inner (outer) cylinder. The extracted mean plume velocity (see text for method) is shown as a 
black line. (b,e,h,k) As above, for every height normalised by the total height of the cylinder 
L. (c,f,i,l) The inverse of the resulting normalised burst velocity. For a = 0 the inner cylinder 
velocity is 1.1 m/s and for a = 0.2 the inner and outer cylinder velocities are 0.91 m/s and 
—0.36 m/s, respectively. The error bars show the width of the test function at 0.95 times the 
peak value. 


typically creating a higher plume velocity than for the a = 0 case. Whereas for the a = 0 
case there was some ambiguity in the direction of the plumes for the azimuthal velocity, 
for a = 0.2 it exactly follows the roll structure. When looking at the radial velocity, the 
inward plumes are less pronounced than the outward ones, likely due to a tendency for 
the Ur plumes to flow outwards, as could be seen for the a = 0 case. 

In addition to the dependence on the axial position and the rotation ratio, in figure 8 
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Figure 8. Turbulent plumes in the radial velocity for varying Ta and a = 0. (a) and (b) as in 
figure 7. (c) The plume velocity normalised by the inner cylinder velocity (uo = 0). Errors are 
of the same order as in hgure 7(c). 


the dependence of the plume velocity on the driving parameter Ta is shown. With 
increasing driving strength, the velocity of the plumes also increases. When the plume 
velocity is normalised by {ui — Uq) it is however nearly constant over two decades of Ta. 

3.5. Logarithmie veloeity profiles 

It has been proposed that the logarithmic temperature profiles in Rayleigh-Benard 
convection and logarithmic velocity profiles in Taylor-Couette flow are triggered by plume 
ejection (Ahlers et al 2014; van der Poel et al. 2015; Ostilla-Monico et al. 20146). These 
Refs, also support the notion that parts of the boundary layer can be turbulent, while 
others are not, and that the transition to the ultimate regime entails plume emission 
from the full extent of the boundary layer. 

We will now further investigate the effect of the plumes on the logarithmic nature of 
the velocity profiles. As we have seen before, the roll structure that appears for counter¬ 
rotation of a = 0.2 provides strong mean radial flows which cause plume emission away 
from either the inner or outer cylinder. 

In figure 9(a) the angular velocity at the inner cylinder is shown. It can be seen that 
there is a significant variation in the shape of the profiles, and that profiles at the axial 
positions where a strong positive radial flow exists (top lines in figure 9(c)) follow the 
shape of a logaritmic law more closely. This effect is even stronger at the outer cylinder 
in figure 9(b), where profiles at an axial position with negative radial flow correspond 
best to a log-law. 

To quantify this phenomenom, the velocity profiles are fitted with the modified (Gross- 
mann et al. 2014) Prandtl-von Karman law of the wall: = 1//^^ ln^+ -b 5, with both 
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Figure 9. Angular velocity profile dependence on mean radial flow for counter-rotation a = 0.2 
and Ta = 4.2-10^ at varying axial positions (data from figure 6(li)). (a) Angular velocity profiles 
near the inner cylinder (f G [0,1/2]). The angular velocity and wall normal distance y~^ are 
dehned in the caption of hgure 4, the colours and line styles of the prohles correspond to those 
in hgure 6. The gray area is the htting domain y^ G [50, 800], chosen to encompass the start of 
the outer region up to mid-gap for both inner and outer cylinder, (b) Angular velocity prohles 
near the outer cylinder (f G [1/2,1]). (c) Normalised radial velocity Ur as in hgure 6(f). (d) The 
von Karman constant nuj found by htting = 1 /Kuj In y~^ B for y~^ G [50, 800] near the inner 

(IC) and outer cylinder (OC). The horizontal axis shows the normalised wall normal velocity, 
equal to Ur and —Ur for the IC and OC respectively, averaged over r G [1/4, 3/4]. The error 
bars have a total length of four times the standard error of the ht. (e) The standard error of the 
htting parameter nuj. 


and B as htting parameters. It should be noted that, as previously, a global torque 
is used to determine the normalisation of the prohles. The local torque can be different 
to the average torque, following the large-scale roll structures. This causes the imperfect 
matching with = y^ and means that the value of the htted parameter can differ 
from the local value. The general shape of the prohle is not affected, however. This shape, 
or specihcally, the deviation from a log prohle, is quantihed by the standard error of the 
htting parameter cr{K^). A low value of (T{i<i^) corresponds to a prohle which closely 
follows a log-law. In hgure 9(d) it can be seen that for a mean wall normal velocity 
away from the cylinder, the value of more closely approaches the known classical 
value of 0.40. Additionaly, the error of the ht is much smaller for positive values of the 
outward wall normal velocity as compared to negative values. We can summarise that for 
regions with strong radial how away from a cylinder surface, plumes will emit from that 
cylinder, which in turn create logarithmic boundary layers. This conhrms results from 
direct numerical simulations (van der Poel et al. 2015; Ostilla-Monico et al. 20146). 
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4. Conclusions 

In this work, we measured velocity profiles, the wind Reynolds number and characteris¬ 
tics of turbulent plumes in Taylor-Couette flow for a radius ratio of 0.5 and Taylor number 
of up to 6.2 • 10^. The flow profiles show Taylor vortices for a Taylor number smaller than 
Ta ^ 10^. The normalised profiles approach the Prandtl-von Karman log-law, although 
at the highest Ta = 6.2 • 10^ the log layer is not yet fully developed. They are in good 
agreement with DNS data from other work (Chouippe et al. 2014; Ostilla-Monico et al 
2014c; Ostilla Monico et al. 2015). Due to the strong curvature of this rj = 0.5 setup, a 
large difference between the azimuthal velocity u~^ and the angular velocity co’+ arises. 
The angular velocity uj~^ resembles a log-law more closely, as suggested by Grossmann 
et al. (2014), based on their Navier-Stokes based theoretical considerations. 

Because of the late onset of the ultimate regime for = 0.5, the measurements with Ta 
up to 6.2 • 10^ are in the classical turbulent regime. For the first time the wind Reynolds 
number has been measured in the classical regime of Taylor-Couette flow and it indeed 
follows the theoretically predicted classical scaling of Rcw oc Ta^^^. 

Moreover, we focused on the interplay between rolls, turbulent plumes and logarithmic 
velocity profiles. At a strong driving of Ta = 4.2 • 10^ no significant coherent structures 
exists for pure inner cylinder rotation, but roll structures appear for counter-rotation at 
a = —ujol^i = 0.2. This behaviour has been observed previously for different values of 
T] and Ta, and shows the important part that rolls play in momentum transfer between 
the two cylinders. 

For inner cylinder rotation only, strong outward plumes are visible in the radial velocity. 
In the azimuthal velocity, the plumes mainly go inwards. The exact mechanism causing 
this difference is yet to be elucidated. For counter-rotation, the roll structures strongly 
influence the direction of the plumes. There is a direct correspondence between the 
direction of the mean radial flow and the direction of the plumes at a certain position. 
The plume velocity in the radial flow profiles increases with Taylor number, and has an 
approximately constant value of one quarter of the inner cylinder velocity for pure inner 
cylinder rotation. 

Lastly, by quantifying the correspondence of the angular velocity profiles to a log-law 
for several axial positions, it was found that in regions with strong radial flow away from 
a cylinder surface, plumes will emit from that cylinder, which in turn create logarithmic 
boundary layers. 

This work confirms predictions about velocity profiles and the scaling of the wind 
Reynolds number and sheds new light on the characteristics and the role of plumes in 
Taylor-Couette flow for a radius ratio r] = 0.5, and hopefully will spark further research 
into these intricate phenomena. 


Acknowledgments 

We thank the technicians of the Department of Aerodynamics and Fluid Mechanics 
at the BTU for their technical support. This work was financially supported by the 
European High-Performance Infrastructures in Turbulence (EuHIT), an European Re¬ 
search Council (ERC) Advanced Grant, and the Simon Stevin Prize of the Technology 
Eoundation STW of The Netherlands. 


REFERENCES 

Ahlers, G., Bodenschatz, E. & He, X. 2014 Logarithmic temperature profiles of turbulent 



Taylor-Couette turbulence at radius ratio 77 = 0.5 


15 


Rayleigh-Benard convection in the classical and ultimate state for a Prandtl number of 
0.8. J. Fluid Mech. 758, 436-467. 

Andereck, C. D., Liu, S. S. & Swinney, H. L. 1986 Flow regimes in a circular Couette 
system with independently rotating cylinders. J. Fluid Mech. 164, 155-183. 

Chouippe, a., Climent, E., Legendre, D. & Gabillet, C. 2014 Numerical simulation of 
bubble dispersion in turbulent Taylor-Couette flow. Phys. Fluids 26 (4), 043304. 

Couette, M. 1890 Etudes sur le frottement des liquides. Ann. Chim. Phys. 21, 433. 

Dong, S. 2007 Direct numerical simulation of turbulent Taylor-Couette flow. J. Fluid Mech. 
587, 373-393. 

Donnelly, R. J. 1991 Taylor-Couette flow: the early days. Phys. Today 44(11), 32-39. 

Dubrulle, B., Daughot, O., Daviaud, E., Longgaretti, P. Y., Righard, D. & Zahn, 
J. P. 2005 Stability and turbulent transport in Taylor-Couette flow from analysis of 
experimental data. Phys. Fluids 17, 095103. 

Egkhardt, B., Grossmann, S. & Lohse, D. 2007 Torque scaling in turbulent Taylor-Couette 
flow between independently rotating cylinders. J. Fluid Mech. 581, 221-250. 

Eardin, M. a., Perge, C. & Taberlet, N. 2014 “The hydrogen atom of fluid dynamics” - 
Introduction to the Taylor-Couette flow for Soft Matter scientists. Soft Matter 10, 3523- 
3535. 

Eenstermagher, P. R., Swinney, Harry L. & Gollub, J. P. 1979 Dynamical instabilities 
and the transition to chaotic Taylor vortex flow. J. Fluid Mech. 94 (01), 103-128. 

VAN Gils, D. P. M., Huisman, S. G., Grossmann, S., Sun, C. & Lohse, D. 2012 Optimal 
Taylor-Couette turbulence. J. Fluid Mech. 706, 118-149. 

Grossmann, S. & Lohse, D. 2000 Scaling in thermal convection: A unifying theory. J. Fluid. 
Mech. 407, 27-56. 

Grossmann, S. & Lohse, D. 2011 Multiple scaling in the ultimate regime of thermal convection. 
Phys. Fluids 23, 045108. 

Grossmann, S., Lohse, D. & Sun, C. 2014 Velocity profiles in strongly turbulent Taylor- 
Couette flow. Phys. Fluids 26, 025114. 

Grossmann, S., Lohse, D. & Sun, C. 2016 High Reynolds number Taylor-Couette turbulence. 
Annu. Rev. Fluid Mech. 48, 53-80. 

He, X., Eunesghilling, D., Nobagh, H., Bodensghatz, E. & Ahlers, G. 2012 Transition to 
the ultimate state of turbulent Rayleigh-Benard convection. Phys. Rev. Lett. 108, 024502. 

VAN Hout, R. & Katz, J. 2011 Measurements of mean flow and turbulence characteristics 
in high-Reynolds number counter-rotating Taylor-Couette flow. Phys. Fluids 23 (10), 
105102. 

Huisman, S. G., van Gils, D. P. M., Grossmann, S., Sun, C. & Lohse, D. 2012 Ultimate 
turbulent Taylor-Couette flow. Phys. Rev. Lett. 108, 024501. 

Huisman, S. G., Sgharnowski, S., Cierpka, C., Kahler, C., Lohse, D. & Sun, C. 2013 
Logarithmic boundary layers in strong Taylor-Couette turbulence. Phys. Rev. Lett. 110, 
264501. 

Huisman, S. G., van der Veen, R. C. A., Sun, C. V Lohse, D. 2014 Multiple states in 
highly turbulent Taylor-Couette flow. Nat. Commun. 5, 3820. 

Ji, H., Burin, M., Sghartman, E. & Goodman, J. 2006 Hydrodynamic turbulence cannot 
transport angular momentum effectively in astrophysical disks. Nature 444, 343-346. 

Kieeer, j. 1953 Sequential minimax search for a maximum. Proc. Amer. Math. Soc. 4, 502-506. 

Kraighnan, R. H. 1962 Turbulent thermal convection at arbritrary Prandtl number. Phys. 
Fluids 5, 1374-1389. 

Lathrop, D. P., Eineberg, Jay & Swinney, H. S. 1992 Turbulent flow between concentric 
rotating cylinders at large Reynolds numbers. Phys. Rev. Lett. 68, 1515-1518. 

Lewis, G. S. & Swinney, H. L. 1999 Velocity structure functions, scaling, and transitions in 
high-Reynolds-number Couette-Taylor flow. Phys. Rev. E 59, 5457-5467. 

Mallogk, a. 1896 Experiments on Eluid Viscosity. Phil. Trans. R. Soc. Lond. A 187, 41-56. 

Marusig, L, Monty, J. P., Hultmark, M. & Smits, A. J. 2013 On the logarithmic region 
in wall turbulence. J. Fluid. Mech. 716. 

Merbold, S., Braugkmann, H. j. V Egbers, C. 2013 Torque measurements and numerical 



16 


R.C.A. van der Veen and others 


determination in differentially rotating wide gap Taylor-Couette flow. Phys. Rev. E 87, 
023014. 

Merbold, S., Fischer, S. & Egbers, C. 2011 Torque scaling in Taylor-Couette flow — an 
experimental investigation. J. Phys.: Conf. Ser. 318 (8), 082017. 

Ostilla-Monico, R., Huisman, S. G., Jannink, T. J. G., van Gils, D. P. M., Verzicco, 
R., Grossmann, S., Sun, G. & Lohse, D. 2014a Optimal Taylor-Gouette flow: radius 
ratio dependence. J. Fluid Mech. 747, 1-29. 

Ostilla-Monico, R., van der Poel, E. P., Verzicco, R., Grossmann, S. V Lohse, D. 
20146 Boundary layer dynamics at the transition between the classical and the ultimate 
regime of Taylor-Gouette flow. Phys. Fluids 26 (1), 015114. 

Ostilla-Monico, R., van der Poel, E. P., Verzicco, R., Grossmann, S. V Lohse, D. 
2014c Exploring the phase diagram of fully turbulent Taylor-Gouette flow. J. Fluid Mech. 
761, 1-26. 

OsTiLLA Monico, R., Verzicco, R., Grossmann, S. V Lohse, D. 2015 The effect of Taylor 
rolls on Taylor-Gouette flow. ArXiv e-prints , arXiv: 1501.03012. 

VAN DER Poel, E. P., Ostilla-Monico, R., Verzicco, R., S., Grossmann V Lohse, D. 
2015 Logarithmic mean temperature profiles and their connection to plume emissions in 
turbulent Rayleigh-Benard convection. Phys. Rev. Lett. Submitted. 

Ravelet, E., Deleos, R. V Westerweel, J. 2010 Influence of global rotation and Reynolds 
number on the large-scale features of a turbulent Taylor-Gouette flow. Phys. Fluids 22 (5), 
055103. 

Taylor, G. I. 1923 Experiments on the motion of solid bodies in rotating fluids. Proc. R. Soc. 
Land. A 104, 213-218. 

Tokgoz, S., Elsinga, G. E., Deleos, R. V Westerweel, J. 2011 Experimental investigation 
of torque scaling and coherent structures in turbulent Taylor-Gouette flow. J. Phys.: Conf. 
Ser. 318 (8), 082018. 

Wendt, E. 1933 Turbulente Stromungen zwischen zwei rotierenden konaxialen Zylindern. 
Ingenieur-Archiv 4, 577-595. 



